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SUMMARY 


Eliminating the corporate feed network in shaped countour beam antennas 
will reduce the expense, weight, and RF loss of the antenna system. One way of 
producing a shaped contour beam without using a feed network is to use a single 
shaped reflector with a single feed element. For a prescribed contour beam and 
feed, an optimization method for designing the reflector shape is given. As a 
design example, a shaped reflector is designed to produce a continental U. S. 
coverage (CONUS) beam. The RF performance of the shaped reflector is then 
verified by physical optics. 


*Work supported by NASA Grant NAG3-419. 



1. Introduction 


In many applications, spacecraft antennas require radiation patterns to be 
shaped such that the pattern contour fits the shape of the desired coverage 
region. A common example is to design the constant gain levels of the pattern 
to follow the projected boundary of the continental United States (CONUS) as 
viewed from synchronous orbit. The shaped contour radiation pattern reduces 
wasted transmitted power by minimizing the illumination of unwanted areas such 
as the ocean in this example. 

The two most popular designs for producing shaped contour radiation pat- 
terns are the array fed parabolic reflector [1], [2] and the direct radiating 
planar array [3] (see Figures la and lb, respectively). If the particular appli 
cation calls for a single fixed shaped beam (no electronic scanning is needed), 
then both approaches generally employ passive beam-forming networks to properly 
weight the array elements. There are, however, several disadvantages to passive 
beam-forming networks. One drawback is that a major part of the cost of the 
antenna is in the construction and tuning of the beam-forming network; construc- 
tion and tuning become more difficult as the frequency is increased. Another 
disadvantage is the RF loss associated with passive beam-forming networks. In 
general, the network loss also becomes larger with increased frequency. 

It is, however, possible to generate a shaped contour radiation pattern 
without using a beam-forming network. This is accomplished by using a single 
feed element and a single shaped reflector (see Figure lc). In this case, the 
far-field radiation pattern is configured to the desired shape by properly 
shaping the reflector surface. Using this design, a relatively inexpensive, low 
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loss shaped beam antenna can be produced with no electrical tuning other than 
alignment of the feed element and reflector. 

The subject of single reflector shaping for arbitrarily shaped contour beams 
has been an active research area since the beginning of the 1970's. One of the 
earliest methods consists of a wavefront synthesis technique [4] . In this 
method, the wavefront of the far-field radiation pattern is assumed to be com- 
posed of two parts. The inner part is a spherical wave limited in angular 
extent by the approximate shape of the desired coverage region. The outer part 
is a ruled surface with the contour of the boundary of the inner part as its 
directrix. The wavefront of the feed element is assumed to be spherical. The 
reflector surface is then completely determined from the incident and reflected 
wavefronts by applying the principles of geometric optics [5], [8], [9]. This 
method yields good results provided design parameters (such as the thickness of 
the outer part of the wavefront) are optimized. The design parameter optimiza- 
tion is necessary because there are usually differences between the desired 
radiation pattern and the actual radiation pattern calculated from the generated 
reflector surface. This difference exemplifies the difficulty of handling the 
geometric optics caustic associated with the outer part of the reflected 
wavefront . 

A later method of reflector shaping replaces the wavefront synthesis with a 
more rigorous aperture phase synthesis [6]. In this method, the aperture phase 
distribution is optimized (for a fixed aperture amplitude distribution) to shape 
the far-field pattern. The reflector surface is then determined from the incident 
and reflected wavefronts in a manner similar to that given by Katagi [4] . 
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Aperture phase synthesis eliminates the caustic problem and allows more pattern 
control (i.e., sidelobe suppression) than wavefront synthesis. However, one 
problem with aperture phase synthesis is that the aperture amplitude distribu- 
tion and the aperture phase distribution cannot be independently specified (as 
in the case of dual shaped reflectors) [7] . The reflector shape that produces 
the optimized aperture phase distribution will not produce the aperture ampli- 
tude distribution assumed in the phase optimization procedure. Jorgensen in [6] 
uses a fixed Gaussian amplitude distribution and assumes the aperture phase has 
the major contribution to beam shaping. This is not a bad approximation, but it 
begins to break down for more complicated radiation pattern contours by causing 
the radiation pattern calculated from the synthesized reflector surface to be 
different from that for the desired radiation pattern. 

To improve on the previous methods, we must find a reflector design proce- 
dure that produces an aperture amplitude and phase distribution pair that satis- 
fies two conditions: 

1) The aperture amplitude and phase distributions produce the desired 
beam shape. 

2) The aperture amplitude and phase distributions can be generated by a 
single shaped reflector. 

In this paper a method will be presented that insures that the synthesized aper- 
ture amplitude and phase distributions satisfy the above two conditions. 
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The synthesis method is an iterative technique that consists of three basic 
procedures. In the first section of this paper, the iterative technique is 
outlined. In the next three sections, each of the three basic procedures is 
described in detail. The last part of this paper applies the iterative tech- 
nique to the design of a CONUS beam. 
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2. Reflector Shaping Technique 


The reflector shaping technique consists of what can be conveniently 
thought of as the repetition of three basic procedures. Consider the reflector- 
feed-near-field aperture plane diagram shown in Figure 2. The first procedure 
is to optimize the phase distribution in the near-field aperture plane for a 
fixed near-field amplitude distribution. This procedure is similar to the 
method of [6] except that the aperture is divided into an array of square ele- 
ments where each element is assumed to have uniform amplitude and phase. A 
least mean square phase optimization routine is then applied to find the best 
phase distribution for the "array" for a given amplitude distribution. The 
second procedure is to calculate the coordinates of the reflector surface that 
create the given phase distribution. As mentioned previously, the reflector 
surface can be completely determined from the phase distribution by applying the 
principles of geometric optics. The third procedure is to calculate the aper- 
ture amplitude distribution from the shaped reflector surface and feed element 
pattern. The new aperture amplitude distribution will be used in the next 
iteration of the phase optimization procedure. 

Figure 3 depicts the iterative technique for determining the reflector shape. 
Since there must be some starting point in any iterative technique, a perfect 
offset parabolic reflector of an arbitrary projected aperture is chosen. This 
reflector produces a certain amplitude distribution that depends on the feed 
pattern and a uniform phase distribution. The first procedure, phase optimiza- 
tion, is then performed. The result is the parabolic reflector amplitude dis- 
tribution combined with an optimized phase distribution that yields a far-field 
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radiation pattern with the desired shape. Next, the second procedure is applied 
and the reflector surface that produces the optimized phase distribution is 
determined. Once the reflector surface is known, the third procedure is performed 
and the new aperture amplitude distribution is obtained. From the aperture 
amplitude and phase information, the far-field radiation pattern of the shaped 
reflector can be calculated by Fourier transformation. If the far-field pattern 
meets the specifications, then the iterative process is ended; if not, then the 
new aperture amplitude distribution is used in the next iteration of the phase 
optimization procedure. The iterative process is then repeated until the beam 
shape meets the specifications, or the beam contour settles to some particular 
shape. The convergence of the beam shape is not guaranteed, but in practice the 
beam contour usually converges. This is due to the fact that the aperture phase 
is the dominant factor in beam shaping; the new aperture amplitude only serves 
to perturb the beam shaping. Consequently, the aperture phase (and thus ampli- 
tude) changes by smaller amounts as the iterative process progresses. 
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3. Phase Optimization Procedure 


To formulate the phase optimization problem, consider the geometry of 
Figure 2. An aperture plane in front of the reflector surface is assumed to be 
located at z = z^. The amplitude and phase of the electric field that is 
reflected from the reflector surface are calculated on a grid of points in the 
aperture plane using geometric optics [5], [9]. The aperture plane field 
distribution is then discretized by dividing it into square patches, where the 
amplitude and phase of each patch are assumed constant over the extent of the 
patch and equal to the value of the reference polarized electric field (i.e., 

Ey) at the center of the patch. This is essentially a two-dimensional pulse 
function expansion of the reference polarized aperture field, and will allow the 
aperture field to be treated as an array of identical elements. 


The normalized reference polarized electric field component for a given 
direction in the far field can be written as 


E(e,*) - l A. E. (6,*)e 
i-1 


jk(x a £u + y ai v) e j6i 


where E^(0,<}>) = the normalized reference polarized electric field component (in 
the far field) due to the i element in the aperture plane 
array. Unity amplitude and zero phase weighting are assumed. 

= relative amplitude of the i th element 
3^ = relative phase of the i ^ element 
N - number of elements in the aperture 

u = sin 9 cos $ 
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v = sin 6 sin $ 

x . = the x coordinate of the i element 
ai 

y a i - the y coordinate of the i fc ^ element 


Note that the normalization is such that the gain in dBi of the i element is 
given by 

Gain.(e,*) = 10 log 10 ( |E.(0 ,+ ) \ 2 ) (2) 


The gain of the array can then be approximated as 


Gain(0 ,<{> ) 


|ECe ,d> ) 

N 2 

l A. 

* 1 1 
1 = 1 


2 


(3) 


Here it is assumed that the power radiated by each individual isolated element 
can be added to obtain the power radiated by the array. 


If (3) is parameterized such that gain in the direction of interest in 
the far field is denoted by 


Gain( 0 


|E(e.,*.)| 2 |e | 2 

.,*.) = jr — = — ■ G. 

J’*J N N j 


(4) 




l A/ 
1 

1=1 


where 


2. = l A.E. . 
J i=i 1 1J 


36i 


E. . = E.(0.,*.)e jk(XaiU j + y aivj) 
i J ,Y J 
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then a cost function that is always positive, can be defined as 


* = 



G. 

J 



( 5 ) 


where 


GDj = the gain desired in the j t ^ 1 direction of interest 
M * number of far-field directions where gain is to be optimized. 


Minimizing ^ is the criterion for the optimization procedure. Note that when ij/ 
is minimized, the mean square difference between the gain desired and the gain 
achieved is minimized. For phase optimization, the set of phases 0^ that mini- 
mize is sought. 

The set of 0^ that minimizes ij/ can be found by applying the method of 
steepest descent. The method takes the following form 


— — 


— 


— 


6 1 (k+l) 


0 1 0c) 


at 

aBj 

Bj » B j(k) 

S 2 (k+1) 


8 2 (k) 


at 

as 2 

B 2 - S 2 (k) 


= 

: 

- V 

• 


0 N (k+l) 


0 N (k) 


at | 

L 38 n 

1 B n = 6 N (k) 


where 0^(k) = 0^ on the k^ iteration of the algorithm of (6). By properly 
choosing p f will be minimized since the gradient (the last term on the right 
of (6)) gives the direction in 0-space that results in the greatest increase in 
(see Figure 4), A method for choosing p that works fairly well in practice 
is to first set p equal to the inverse of the magnitude of the gradient in (6), 
then reduce p by an additional 0.9 times every time the cost function increases 
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The gradient can be easily calculated since ip is an explicit function of 


The result is 


*- - I 2[G - GDj ] -g-L- [j* iSi t * - jAjEy* e' j6i I .] 

1 J“1 J J V A 2 


8 ^ = 
86 


(7) 


y a. 

*i 1 

i=*i 


Note that the gradient is also a function of the particular set of A^, i=l, N 
used in the procedure. 

By starting with initial phase and amplitude distributions, 6^(0) and A^(0) 
i-l,N, a new phase distribution, 6^(k) i-l,N, (A^O) i*l,N remain unchanged) is 
obtained that minimizes the mean square difference between the gain desired and 
the gain achieved at the various far-field directions. By choosing the desired 
gains and their corresponding directions to cover the area of interest (CONUS, 
for example), a shaped contour radiation pattern can be produced. Note that 
sidelobes can be suppressed in selected directions by choosing the gain desired 
to be small for those directions. 

The choice of the least mean square optimization criterion (5), (6) is 
arbitrary. The beam shaping method could also employ a minimax algorithm as 
long as phase is the only optimized variable. 
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4. Calculation of the Reflector Surface 


Once the optimized aperture phase distribution is known, a reflector sur- 
face can be completely determined from this distribution by applying the prin- 
ciples of geometric optics. Consider the geometry of Figure 5. If one reflec- 
tion point. A, on the reflector surface is assumed to be known, then the length 
of what will be referred to as the reference path, d, is known and given by 

d = | H | + | Xa s | (8) 

where F denotes the feed position at (0, 0, f) 

A denotes a point in the aperture plane. 

Note that the incident ray t’A produces the reflected -ray %A , and the total 

a 

length along this ray path is given by the reference path length. For any other 
point P on the reflector surface, the length of the ray path d' is given by 

d' = |£p|+|?p| = [x 2 + y 2 + (z-f ) 2 ] 1/2 + [(x-x ) 2 + (y-y ) 2 + (z-z ) 2 ] 1/2 (9) 

3 L 3 3 3 

where P denotes a point in the aperture plane (x , y , z ) 

a a a a 

P denotes a point on the reflector (x, y, z) . 

In this case the incident ray £p produces the reflected ray ?P . If it is 
assumed that the phase in the aperture plane is continuous and that the phases at 
points A and P in the aperture plane are given by 0(A) and 0(P ), respectively, 
then the following relation holds 

d' = — j£- [e(P a ) - e(A a )] + d do) 
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where 


k 


2tt 

X 


Equation (10) relates the length difference between the two ray paths to the 
difference in their corresponding phase values in the aperture plane. The 
right-hand side of (10) is known and the left-hand side contains the unknowns x, 
y, and z, the coordinates of a point on the reflector surface which remains to 
be determined. Two more independent equations are needed to solve for the three 
unknowns. These two equations are obtained from the expression for the line 
passing through the points P and P . In the geometric optics limit, the ray 
$p is normal to the equiphase surface and parallel to the wave vector £ at the 

3 

point P . More precisely 
a 


?P 




-T 

k 

k 

a 7 (k x+k y + k 2) 
k ^ x y 7 z * 

= i ve (x, y, z) 


P 

P 


a 

a 


(id 


where 


2tt 

A 


1*1 - * - 

Therefore, the two equations for the line passing through ?P are given by 


X-X Z-Z 

a a 


m 


m 


( 12 ) 


y-y a 


m m 

y 2 


(13) 


where 


m = i — 
x k 


I i! 

k 9x 


( 14a) 


k 

- -t L 

y k 


_ J_ 39_ 
k 3y 


(14b) 
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m z - [i -m x 2 - m y 2 ] 1/2 (14c) 

Since the phase distribution in the aperture plane is known, the partial deriva- 
tives of the phase with respect to x and y can be approximated by difference 
expressions. The only unknowns in Equations (12) and (13) are x, y, and z. 
Therefore, Equations (10), (12) and (13) can then be solved for the three 
unknowns (x, y, z) yielding a point on the reflector surface. The solution of 
the three equations is given by 


z = 


2 2 2 2 2 
Q (x a * y a ) ~ x a ~ y a ~ *a * 2z a f ~ f 

2(cx + yy + z - f) - 2Q(x , y )(o 2 + y 2 + 1) 

d a a a a 


1/2 " Z a 


(15a) 


y = y (z ~ z a ) + y a 


(15b) 


x = cr (z - z ) + x 
where * 

Q<«,. y*> l«< A a > - e<V) * d 



(15c) 


m 


x 

a = — 
m z 

At the beginning of thi3 derivation two variables were assumed to be known, 
0(A a ). These variables are chosen arbitrarily as 


d, 


d = f + z (16a) 

a 

6(A ) = 0 (16b) 

a 

where f is the focal length of the initial parabolic reflector that is used to 

begin the iterative surface shaping technique and z = 0 is the equation of the 
' a 
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aperture plane. Note also that the solutions of (10), (12) and (13) automati- 
cally satisfy Snell’s law for smooth phase distributions [8] • 

It should be mentioned that if there are discontinuities in the aperture 
phase distribution, the corresponding reflector surface points will not satisfy 
Snell’s law. Since, for every point in the aperture plane, a reflection point 
on the reflector surface is determined, a normal for every reflector point can 
be determined from neighboring reflection points. In this way, Snell’s law can 
be tested. It will be shown that discontinuities in the aperture phase will 
produce discontinuities in the reflector surface. 
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5. Calculation of the Aperture Plane Amplitude Distribution 


After the surface has been determined at a discrete set of points, the 
aperture amplitude distribution can be calculated. All points on the reflector 
surface can be tested to see if Snell's law is satisfied. For those points that 
satisfy Snell's law, the amplitude of the exit ray can be calculated in the 
aperture plane by the projected area ratio of a ray tube (see Figure 6). It is 
well known [5], [9] that 




/Area 
/ Area 


1 

2 


(17) 


where |E r | = magnitude of the reflected field at the reflector surface 

|E | = magnitude of the reflected field in the aperture plane 
Areaj = projected area of the ray tube at the reflector surface 
Area 2 = projected area of the ray tube in the aperture plane. 


A small ray tube can be formed by selecting three rays in close proximity. 

Consider Figure 6; the three incident rays are $P(I,J), ?P(I+1,J) and 

fp(I,J+l). The corresponding reflected rays are P( I , J)£ (I,J), P( 1+1 , J)P (I+1,J) 

a a 

and P(I,J+1)£ (I,J+1). As the aperture points and the reflection points are 
brought closer together, (17) can be used to approximate the aperture field at 
point P (I,J). The approximation becomes exact as the areas go to zero. The 
projected areas are given by (see Figure 7) 

Area^ = y h • side 1 * to • n (18) 

where h - side 2 * / l-(cos a)^ 

= (side 3) 2 - (side l) 2 - (side 2) 2 
C0S ° 2(side 1) * (side 2) 
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side 1 = | P(I,J) - P(I+1,J) | 


The field 
therefore , 
(17). 


side 2 » | P(I,J) - P(I,J+1) | 
side 3 - | P(I+1,J) - P(I,J+1) | 

to * normalized incident ray direction 
n = normal to the surface at P(I,J) 

Area 2 58 y * cell * fc • z 

cell = area of square patch in aperture 

^r * normalized reflected ray direction 
z = normal to aperture plane 

quantity |E^ | in (17) is known from the given primary feed pattern 
the amplitude distribution in the aperture can be calculated from 


(19) 
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6. Design of a CONUS Beam 


The iterative technique outlined in the previous sections was used to 
design a CONUS shaped contour radiation pattern. To begin the design procedure, 
73 directions in the far field were chosen to cover the projected CONUS map as 
viewed from a particular synchronous orbit position as shown in Figure 8. Next, 
an initial offset parabolic reflector and feed element were selected for the 
starting point of the iterative process. This reflector and feed are shown in 
Figure 9, and the far-field pattern is shown in Figure 10. The amplitude and 
phase of the reference polarized electric field component were then calculated 
over the aperture plane on a 47 by 47 point square grid. The resulting 2209 
square patches (elements) in the aperture plane were each 0.543 X on a side. 

The normalized reference polarized electric field component due to each patch 
was calculated in each of the 73 far-field directions. The field component due 
to each of the 2209 aperture elements in each of the 73 far-field directions 
forms the data base for the phase optimization procedure. 

The three basic procedures - phase optimization, reflector surface calcula- 
tion, and aperture amplitude distribution calculation - were performed using the 
initial aperture amplitude and phase distributions for the perfect parabolic 
reflector. From the resultant new aperture amplitude and phase distributions 
(and corresponding reflector surface), the far-field pattern shown in Figure 11 
was obtained. The beam shape is beginning to take on the CONUS form, but the 
process needs to be continued to improve the beam shape. 

The three basic procedures were then repeated for a second iteration with 
the aperture amplitude and phase that produced the far-field pattern shown in 
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Figure 11 as input. The resultant aperture amplitude and phase after the second 
iteration produced the far-field pattern shown in Figure 12. The procedures were 
again repeated for a third iteration, and the far-field pattern of Figure 13 was 
produced from the resultant shaped reflector surface (and corresponding aperture 
amplitude and phase). After the fourth iteration, the far-field pattern of 
Figure 14 was obtained. The far-field pattern did not change much when going 
from the third to the fourth iteration so the iterative process was ended. The 
final aperture amplitude distribution and aperture phase distribution are shown 
in Figures 15 and 16, respectively. The final shaped reflector surface is shown 
in Figures 17 and 18. The cross-polarization pattern after the fourth iteration 
is shown in Figure 19. (The feed model had no cross-polarized component.) The 
gain over the coverage region was greater than 28.0 dBi with the polarization 
isolation less than -26 dB. 

Note that the difference between the shaped surface and the parabolic sur- 
face is almost saddle shaped with slight discontinuities in the surface in the 
steep portions of the saddle shape. Notice similar discontinuities in the phase 
distribution. These discontinuities were detected by testing to see if Snell's 
law was satisfied. At the discontinuity points (edge points) the reflected 
fields were neglected when calculating the far-field pattern. The electrical 
performance of the shaped reflector was checked by applying physical optics. 
Figure 20 shows the physical optics far-field radiation pattern for the inter- 
polated shaped reflector surface of Figure 18. Note there is very little dif- 
ference between the physical optics pattern and the geometric optics pattern of 
Figure 14. 
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7. Conclusion 


A method has been presented for synthesizing a single reflecting surface 
that accurately produces a shaped contour radiation pattern when the radiation 
contour shape and feed pattern are specified. The method is quite general in 
that it will take into account discontinuities in the reflecting surface. 
Results were generated for a CONUS beam design, and very good beam shaping 
resulted. The electrical performance of the shaped reflector surface was also 
verified by physical optics. 
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Beam Forming Network 


FIGURE 1b. Direct Radiating Planar Array 


23 







FIGURE 2. Geometry of the Reflector Shaping Problem 
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FIGURE 3. Iterative Reflector Shaping Technique 







figure 7. Cross Sections of the Ray Tubes at the Reflector 
and In the Aperture Plane 
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FIGURE 8. Farfleld Points Over which Gain Is Optimized 
(Desired Gain > 28 dBi for all Points) 
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FIGURE 9. Geometry of the Initial Parabolic Reflector 
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FIGURE 10. Far Field Pattern for the Initial Parabolic Reflector 
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FIGURE 11. Far Field Pattern for the Shaped Reflector after the First Iteration 
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FIGURE 12. Far Field Pattern for the Shaped Reflector after the Second Iteration 




FIGURE 13. Far Field Pattern for the Shaped Reflector after the Third Iteration 
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FIGURE 14. Far Field Pattern for the Shaped reflector after the Fourth 
Iteration. The antenna gain is greater than 28.0 dBi over the 
coverage area. 



FIGURE 15. Aperture Amplitude Distribution after the 
Fourth Iteration 
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FIGURE 16. Aperture Phase Distribution after the 
Fourth Iteration 
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FIGURE 17. A Plot of the Difference (in z Values) Between 
the Shaped Reflector Surface and the Initial 
Parabolic Reflector Surface. The shaped surface 
data was interpolated on a uniform x-y grid from the 
randomly spaced data obtained from the numerical results. 

The peak to peak variation is about 0.5 X. 
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FIGURE 20. Physical Optics Farfield Pattern for the Shaped Reflector. 

The pattern was calculated using physical optics and the 
interpolated reflector surface after the fourth iteration. 
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